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Abstract 

According to band structure calculations, the Fermi surface of the quasi-two dimensional metal 
0-(ET)4ZnBr4(C6H4Cl2) illustrates the linear chain of coupled orbits model. Accordingly, de Haas- 
van Alphen oscillations spectra recorded in pulsed magnetic field of up to 55 T evidence many 
Fourier components, the frequency of which are linear combinations of the frequencies relevant to 
the closed a and the magnetic breakdown /3 orbits. The field and temperature dependence of these 
components’ amplitude are quantitatively accounted for by analytic calculations including, beyond 
the Lifshitz-Kosevich formula, second order terms in damping factors due to the oscillation of the 
chemical potential as the magnetic field varies. Whereas these second order terms are negligible 
for the orbits a, (3 and 2/3 — a, they are solely responsible for the ’forbidden orbit’ /3 — a and its 
harmonic and have a significant influence on Fourier components such as 2a and /3 -|- a, yielding 
strongly non-Lifshitz-Kosevich behaviour in the latter case. 

Keywords: De Haas-Van Alphen oscillations, high magnetic fields, two-dimensional organic met¬ 
als. 
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I. INTRODUCTION 


Many quasi-two-dimensional (q- 2 D) metals, in particular charge transfer salts based on 
the bis-ethylenedithio-tetrathiafulvalene (ET) molecule, illustrate the textbook Fermi sur¬ 
face (FS) proposed by Pippard in the early sixties. This model FS, an example of which is 
provided in the inset of Fig. [H was intended to compute the de Haas-van Alphen (dHvA) 
oscillation spectrum of the linear chain of orbits coupled by magnetic breakdown (MB)^!^. 
In line with the coupled orbits network model of Falicov-Stachowiak->^, relevant dHvA oscil¬ 
lations spectra involve linear combinations of frequencies linked to the a and MB-induced f 3 
orbits^'-. However, it is now well established that the held and temperature dependence of 
many of these Fourier components cannot be accounted for by this model due to oscillation 
of the chemical potential in magnetic heldl^^— . Analytic tools, given in the Appendix, have 
been provided in order to quantitatively account for the held- and temperature-dependent 
amplitudes of the various Fourier components observed^^— . Briehy, in addition to a hrst 
order term corresponding to the Lifshitz-Kosevich (LK) model^, second order terms due to 
oscillation of the chemical potential must be taken into account. Nevertheless, their rela¬ 
tive importance strongly depends on the involved parameter values, in particular the Lande 
factors. As an example, provided spin damping factors relevant to basic orbits are not too 
small, i. e. gljUa and (where and ma(i3) are the ehective Lande factor and ehec- 

tive mass, respectively, of the a{( 3 ) orbit) are not close to odd integers, these second order 
terms have a negligible contribution to the Fourier amplitude and Ap, respectively. In 
contrast, the amplitudes Ap(^p_a) of the Fourier components with frequencies p(Fs — Fk), 
which are commonly referred to as ’forbidden orbits’ since they do not correspond to MB 
orbits, are only governed by second order terms. For completeness. It should be noticed 
that, in the case of magnetoresistance oscillations, components such as /3 — a or /3 — 2 a 
correspond to quantum interference paths^^ which are liable to enter the Shubnikov-de Haas 
(SdH) spectrai^i. 

Up to now, these calculations have only been implemented to account for the data of the 
strongly two-dimensional compound 0-(ET)4CoBr4(C6H4Cl2)^'^, referred in the following 
to as the Co-compound. For this compound, the held and temperature dependence of 
the second harmonic amplitude of the a orbit (A2Q,), which signihcantly differs from the 
predictions of the LK model, and the ’forbidden orbit’ amplitude are quantitatively 


2 


accounted for by the calculations. Nevertheless, data analysis for other compounds, with 
different FS parameters, are needed to further check the model. In addition, depending on 
the value of the involved FS parameters (in particular effective masses and Lande factors), 
strongly non-monotonic held and temperature dependence is liable to be observed in few 
cases^il^. Actually, such a feature has never been reported yet. 

The aim of this article is to report on quantum oscillations spectra of 0-(ET)4ZnBr4(C6H4Cl2), 
referred to as the Zn-compound in the following. This compound belongs to the same family 
as the Co-compound, namely ^-(ET) 4 MBr 4 (CfiH 4 Ch), where M is a metal such as Co, Zn, 

Hg, Cd (for a review, see Ref. l22(). Strikingly, the crystal structure of these compounds 
involves one conducting and one insulating ET plane, with different atomic arrangement, 
insuring a strong two-dimensionality. More extended data than for the previously reported 
Co-compound, i.e. held and temperature dependence of Fourier amplitude relevant to sev¬ 
eral frequency combinations, are derived, allowing a more extensive check of the formulas 
reported in the Appendix. In particular, it is demonstrated that strongly non-monotonic 
temperature dependence of the Fourier component with frequency corresponding to the MB 
orbit fi+a is observed. 


II. EXPERIMENTAL 


Crystals were synthesized by electrocrystallization technique as reported in Ref. 
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The FS topology was obtained through extended Hiickel type tight-binding band struc 


ture calculations^i^ as reported in Refs. 


25 


and 


26 


These calculations were based on X-ray 


dihraction data collected at 100 K and 180 K at the IPCP-Chernogolovka and the LCC- 
Toulouse, respectively. 

Six crystals denoted hereafter as crystal #1 to #6, respectively, were studied in pulsed 
magnetic helds of up to 55 T with a pulse decay duration of 0.32 s. DHvA oscillations were 
measured through magnetic torque measurements of crystals #1 to 7 (^ 4 , with approximate 
dimensions 0.1 x 0.1 x 0.04 mm^, stuck on a microcantilever. Variations of the micro¬ 
cantilever piezoresistance were measured at liquid helium temperatures with a Wheatstone 
bridge with an ac excitation at a frequency of 63 kHz. The angle between the normal to the 
conducting plane and the magnetic held direction was 0 = 11°, 8 ° and 9° for crystals 7 (^ 1 , 7)^2 
and 7 )^ 3 , respectively, while 9 was varied from 15° to 71° thanks to a rotating sample holder 
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FIG. 1. (color on line) (a) Oscillatory part of the high field range of TDO frequency (crystal ^6), 
4-point interlayer magnetoresistance (crystal #5) and torque (crystal ^3) data at 2 K and (b) 
corresponding Fourier analysis (Fourier spectra are shifted down from each other by a constant 
amount for clarity). Thin lines in (b) are marks calculated with Fa =0.93 kT and Fa/Fjs = 0.205. 
The inset displays the Fermi surface of the conducting layer at 100 K. Green rectangle depicts the 
first Brillouin zone where S = {—a*/2, h* Y = (0, 6*/2), M = (a*/2, b*/2), X = (o*/2, 0) and 
F = (0, 0). The basic orbits a and [3 are marked by the blue and red lines, respectively. 
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for crystal ^^4. SdH oscillations were measured through 4-point interlayer magnetoresistance 
(crystal and contactless tunnel diode oscillator (TDO)-based methodl^'^ (crystal #6). 


III. RESULTS AND DISCUSSION 

In the next section fllll Ap the oscillatory data are examined at the light of band structure 
calculations. Section IIII Bl reports on the field and temperature dependence of the basic 
Fourier components, linked to the a and {3 orbits, observed in dHvA and SdH spectra while 
dHvA frequency combinations are considered in Section fill Cl 


A. Band structure calculations and oscillatory spectrum 


Crystalline and electronic band structures of the Zn-compound are very similar to those 
of the Co-compound. Briefly, two different cation layers, labeled A and B, respectively, in 
Refs.— with different atomic arrangements, are observed within the unit cell. According 
to band structure calculations, layer A with a-type packing is insulating while layer B with 
6^-ty pe packing is conducting (for details regarding atomic packing in organic metals, see 


Ref. 


28l ). As reported in Fig. [T], the FS topology relevant to layer B illustrates the Pippard’s 


model, observed in many organic conductors based on the FT molecule. Namely, it is 
composed of one q-2D closed tube and two q-lD sheets separated by a gap. In magnetic 
fields, the closed tube yields the a orbit while, thanks to MB, the jS orbit with an area 
equal to that of the first Brillouin zone (FBZ) is observed. The area of the a orbit is 17.0 
% and 18.2% of the FBZ area at 180 K and 100 K, respectively. It can be remarked that 
the FS of Fig. [1] differs from that of other 0-phase salts. In these latter salts the anions 
impose a periodicity along the b direction which is different from that observed in the Zn- 
and Co-compounds, yielding different FS topology^. 

Fig. [IKa) displays oscillatory parts of the magnetic torque, 4-point longitudinal magne¬ 
toresistance and TDO data at 2 K. It can be remarked first that, while TDO and 4-point 
magnetoresistance data are in phase, magnetic torque data are phase-shifted by 7r/2. This 
feature indicates that, while magnetic torque yields dHvA oscillations, both TDO and 4- 
point magnetoresistance yield SdH oscillations, in agreement with previous statement8^>^. 
Corresponding Fourier analysis are displayed in Fig. [I](b). The two main frequencies, 7% = 
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FIG. 2. (color on line) Temperature dependence (mass plots) of the magnetic torque amplitudes 
(a) and (b) Solid lines in (a) and (b) are best fits of Eos. lA.ll and IA.31 respectively, to the 
data. They are obtained with ma= 1.85, mp= 3.4, Bq = 26 T and Td = 0.8 K. The considered 
magnetic field values are evenly spaced in 1/B in the explored field range, the boundary of which 
are indicated in the figures. 

0.930(2) kT and F/j = 4.534(7) kT, correspond to the a and (3 orbits, respectively, hence 
the a orbit area amounts to 20.5 % of the FBZ area. This value is in good agreement with 
the above reported band structure calculations which are based on X-ray diffraction data 
measured at higher temperature, owing to the increase of the closed tube area relatively 
to that of the FBZ, as the temperature decreases from 180 K to 100 K. Strikingly, an un¬ 
precedented number of frequencies is observed, in particular in the case of the TDO data, 
accounting for the strong non-sinusoidal oscillatory part of these data. These frequencies, 
labeled Frf in the following, are linear combinations of Fa and F^. Frequency as high as 17.2 
kT, corresponding to r; = 4/3 — a, is observed in the TDO spectrum of Fig. [TJ 

B. Basic Fourier components amplitude 

Let us consider hrst magnetic torque data, the oscillation spectra of which involve Fourier 
components with frequencies and amplitudes A/). Since we are dealing with a 2D metal, 
these amplitudes are related to dHvA oscillations amplitudes Arj as oc A/)/Btan(6') where 
6 is the angle between the magnetic held direction and the normal to the conducting plane. 
According to the LK formula, ln(A/)/T) is predicted to vary linearly with the temperature 
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FIG. 3. (color on line) Angle dependence of the magnetic torque amplitudes and Solid lines 
are best fits of Eqs. lA.ll and Eq. IA.3I to the data for and respectively. They are obtained 
with the same effective mass and MB field as in Fig. [2] and 9a — 9’^ — 1-85. 
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FIG. 4. (color on line) Field dependence of the effective mass value of a and /3 orbits derived from 
magnetic torque, 4-point magnetoresistance and TDO data. Horizontal lines mark the effective 
mass values deduced from magnetic torque data, yielding dHvA oscillations, nia = 1-85 and = 

3.4. 
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at a given magnetic field value (mass plot) at high T/B ratio. Since crystals ^1, 41^2 and 
7)^3 yield same results within the error bars reported below, we focus in the following on 
the data relevant to crystal 7 )^ 3 , with the lowest Dingle temperature. Data for rj = a and 
(3 are reported in Fig. [2J They can be analyzed through Eqs. lA.ll and IA.31 respectively. 
However, no fewer than seven parameters enter these equations: effective masses rriadB), 
Dingle temperatures ToaiP), effective Lande factors ga{i 3 ) and MB held Bq. Nevertheless, 
as observed in the case of the Co-compoundl^ and discussed in Ref. Il8|, the second order 


terms of Eqs. lA.ll and IA.3I are negligibly small compared to their leading terms, provided 
the spin damping factors Ra{i 3 ) are far enough from spin-zeroes. As a result, the LK model 
applies and the spin damping factors act as held- and temperature-independent prefactors. 
However, in addition to the ehective masses, the MB held Bq and the Dingle temperatures 
ToaiiB) govern the held dependence of Aa(/ 3 )- As a result, each of the two equations lA.ll 
and IA.3l still involve 3 parameters yielding large uncertainties. For this reason, it is assumed 
in the following that Toa = Within this assumption, data yield ma= 1.85(10), m/ 3 = 
3.40(15) and Bq = 26(3) T for all the three studied crystals. The Dingle temperature, which 
is the only crystal-dependent parameter is T^i = 0.9(1) K, Td 2 = 1-1(1) K and = 
0.8(1) K for crystal 7 )^ 1 , #2 and #3, respectively. Ehective Lande factors which are 

the remaining parameters to be determined, are obtained through the angle dependence of 
Aa(/3). Solid lines in Fig. [3] are the best hts of Eqs. lA.ll and IA.3I to the data relevant to 
crystal #4, yielding g*^ = g* = 1.85(10). 

In short, the ehective masses and Dingle temperatures of the Zn-compound are close 
to the data obtained for the Co-compound whereas the MB held of the latter is higher—. 


Q 


the second 


Owing to the ehective Lande factors values, which were only estimated in Ref. 
order terms of Eqs. lA.ll and IA.3I are negligible in the held and temperature range explored, 
indicating that the LK model, i.e. the hrst order term of Eqs. lA.ll and IA.3I satisfactorily 
accounts for the basic orbits a and (3, respectively. 

A hallmark of the validity of the LK formula is the held-independency of the ehective 
mass derived through this formula from the temperature dependence of the amplitude, as 
it can be observed in Fig. |4] for the dHvA data. In contrast, an apparent strong increase 
of the ehective mass is observed in the case of TDO and, to a less extent, of the 4-point 
interlayer magnetoresistance data. This behaviour can be ascribed to the failure of the 
LK formula for SdH oscillations relevant to basic orbits of q-2D metals^^— . This feature. 





























FIG. 5. (color on line) Mass plot of magnetic torque amplitudes for different values of the magnetic 
field of (a) /3 — a, (b) 2(/3 — a), (c) 2a and (d) 2/3 — a. Solid lines in (a), (b), (c) and (d) are 
calculated with Eos. IA.51 IA.61 and IA.81 respectively. They are obtained with the same set of 
parameters as in Figs. El El Dashed lines in (c) and (d) are obtained with the Lifshitz-Kosevich 
formula. 

which is beyond the scope of the present stndy focused on dHvA spectra, requires specific 
calculations of the conductivity^^— taking into account the multiband nature of the FS. As 
for the smaller discrepancy observed for 4-point magnetoresistance compared to contactless 
TDO measurements, it must be considered that interlayer resistance {Rzz) and in-plane 
resistance {Rxx) which are governed by different matrix elements^, are measured in the 
former and latter case, respectively. Besides, the electrical contacts on the crystal, in the 
case of 4-point magnetoresistance, connect the quasi-particles to a non-quantized reservoir, 
liable to induce a damping of the chemical potential oscillatiort^. 


C. Frequency combinations 

Since the parameters entering Eqs. lA.ll to IA.8I are determined from the analysis of the 
data relevant to the a and /3 orbits, the field and temperature dependence of all the other 
Fourier components amplitude should be accounted for by these parameters. As examples. 
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FIG. 6 . (color on line) Fourier analysis in the frequency range around F/ 3 + 0 , at the mean magnetic 
field value S = 47 T. In the temperature range (a) 2.25-2.56 K and (b) 2.66-3.48 K, the Fourier 
amplitude decreases and increases, respectively, as the temperature increases, (c) Mass plot for 
different values of the magnetic field. Solid lines are calculated with Eq. IA.71 They are obtained 
with the same set of parameters as in Figs. [21 [31 [SJ Dotted lines are the contributions of the first 
order term, corresponding to the Lifshitz-Kosevich formula. 
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the behaviour of few of these amplitudes is considered in Fig. [S] a very good agreement is 
indeed observed. 


Let us examine these data in more details. First, since the f3 — a amplitude is dominated 
by the product Ra,iRp,i (see Eq. IA.5j) its temperature dependence displays a slightly negative 
curvature. Analyzed through the LK formula, the data would yield an apparent effective 
mass close to that oi (3 + a, actually about 0 . 8 (m /3 + nia). This result is in agreement 

with both experimental data relevant to k-ET 2 Cu(NCS )2 for which ~ 0.85(m^ + mQ,)^, 

and in line with numerical simulations^^. 

As already observed^'^'^, 2a is not accounted for by the LK formula as well. Indeed, 
according to Eq. IA.21 the second order term which is of the same order of magnitude as the 
LK damping factor 2 is dominated by R^ ^ accounting for a non-LK behaviour. 

Oppositely, the amplitude of the 2(3 — a component is very close to the prediction of the 
LK model since the second order term, dominated by the product Ra,iRi 3,2 is very small 
compared to the LK term which is proportional to 

Finally, the Fourier component f3 + a is considered in Fig. [HI While below about 2.6 K, the 
amplitude decreases as the temperature decreases, it increases in the range 2.6 K ~ 3.5 K, in 
strong discrepancy with the LK formula. This behaviour is quantitatively well accounted for 
by Eq. lA.71 which evidences a dip in the temperature dependence of the amplitude. Indeed, 
the second order term of Eq. IA.7l is dominated by the product Ra,iRfi,i which contributes to 
the amplitude with an opposite sign to the first order Lifshitz-Kosevich term proportional 
to i?a+/ 3 ,i- These two factors cancel each other at a given held and temperature value (e.g. 
2.9 K at 50 T in the present case), depending on the spin damping factors value, hence on 
the respective values of the products and g^mp. It can be remarked that this feature 

is not observed in the Co-compound^^. Indeed, owing to slightly different effective masses 
and effective Lande factors, the dip in the (3 3- a amplitude would be observed around 9 K, 
i.e. beyond the temperature range in which oscillations can be observed^. 


IV. SUMMARY AND CONCLUSION 

Band structure calculations relevant to the quasi-two dimensional metal 6 '-(ET) 4 ZnBr 4 (C 6 H 4 Cl 2 ) 
indicate that this compound illustrates the linear chain of coupled orbits model proposed 
by Pippardii^ (see Fig. [T]) as it is the case for many organic conductors based on the ET 
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molecule. In line with this statement, quantum oscillations spectra evidence many Fourier 
components, the frequency of which are linear combinations of the frequencies relevant to 
the closed a and the magnetic breakdown (3 orbits. The held and temperature dependence 
of the de Haas-van Alphen amplitude of these components is quantitatively accounted for by 
the analytic calculations reported in the Appendix. Beyond the Lifshitz-Kosevich formula, 
they include second order terms arising from the chemical potential oscillations. These 
second order terms have negligible contributions to the amplitude of the basic a and jS 
components allowing the determination of the various physical parameters entering the data 
(effective masses, magnetic breakdown helds, etc.). They have also a minor contribution on 
the magnetic breakdown orbit 2(3 — a. Oppositely, they have signihcant contribution to 2a 
and (3 + a. Although this latter component physically corresponds to a magnetic breakdown 
orbit, its temperature dependence evidences a strong dip due to the cancelation of the hrst 
and second order terms. Finally, the ’forbidden frequency’ (3 — a and its harmonic 2/3 — 2a, 
which are due to the oscillation of the chemical potential, are accordingly accounted for by 
second order terms, only. 
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Appendix: Analytical expressions of Fourier amplitudes 

In this appendix, are recalled the analytical equations for de Haas-van Alphen amplitudes 
Apri with frequencies pF^ given in Refs.— . They are relevant to two-dimensional FS 
illustrating the Pippard model in which the component is a linear combination of the a 
and (3 orbits and p is the harmonic order (see insert of Fig. [1]). 


A _ T? 

/ia — - ^a,l - 

nrua limp 


-Ra,lRa,2 + -Ra,2Ra,3 + ‘^Rl3,lRa+l3,l + ” -^/3,2-R2^-o,l( A-. 1 ) 
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^/3 D 

-«/3,l- 


-^ 0,1 ~ -^Ra,lRa,3 — Ra,2Ra+^,2 


TrrrijS 

Fp 


TTmp 


-Ri3,iRi3,2 + -Ri3,2Ri3,3 + ‘^Ra,lRa+l3,l + ‘^Rl3,lR2/3,l 


(A.2) 

(A.3) 


[-R/3,2 + 2-R2/3,i] + 




27rm/3 ■ ' ' ' ' vrm^ 

2-Rq,1-R2^-o, 1 — -R/3,2-R2/3,2 “ -R/3,4-R2/3,l] 
Fp-a 


2 2 1 

-^^,1 ~ -^Rp,lRp,3 ~ --R/3,2-R/3,4 — -R«,2-Ra+/3,2 


A 


2(/3-a) 

A^-)-q. 

^2/3-a 


Tirnp 
F2l3-2a 
TTmp 


[Ra,lRl3,l + -Ro,2-R«+/3,1 + RjifiRa+ji,! + -R/3,l-R2/3-a,l] 


2-Ra,2-R2/3,l + 2-R2/3-a,2-R2/3,l + 2 -Rq,,1-R2/3-o, 1 + ^-^«:4-Ro+^,2 


1 1 

-Ra,2Rl3,2 + —R0,2R2l3-a,2 


‘^Fp+a jj Ff^j^a. 

-K/3+q,1 + 


7r?Ti^_i_Q, 

F2p-a 

7rm2p-a 


■R2P-a,l — 


irmp L 
-^2/3-0 


Ra,lRp,l — 2-Ra+/3,2-Ro+/3,l — g-R/3,3-R2^-o,l 


Timp 


-Ra,lRp,2 + -Ra,3Ra+P,2 


(A.4) 

(A.5) 

(A.6) 

(A.7) 

(A.8) 


Damping factors are given by the LK and coupled orbits network models^ as Rri^p{B,T) 
= R'^p{B,T)R^p{B)R^^^{B)R^ .p--^ where the temperature, Dingle, MB and spin damping 
factors are expressed as R^^ = pX^ sinh“^(pA'j,), 7?^^ = exp(—cos 6*)“^), = 

(ipo)”H<?o)”A = cos{Tig*mjj/2cos9), respectively. The held-and temperature-dependent 

variable (X^) and the constant (uq) are expressed as = uorripT/ {B cos 6) and uq = 
27T‘^kBTne{eh)~^ = 14.694 T/K. The tunneling (po) and reflection (go) probabilities are given 
by Po = e“^°AScos6» _|_ g ,2 _ Td is the Dingle temperature dehned by Td = 

h{27rkBT)~^, where is the scattering rate, Bq is the MB held, and g* are the effective 
masses and effective Lande factor, respectively. It can be noticed that the terms of hrst 
order in damping factors correspond to the LK model. 
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